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q ", ABSTRACT 

o 

^ \ We present Infrared Space Observatory spectroscopy of sites in the supernova 

remnants W 28, W 44, and 3C 391, where blast waves are impacting molecular clouds. 
The complete wavelength range from 42 to 188 ^m was observed with the Long 
(--. Wavelength Spectrometer, as well as narrow ranges centered on 4.695, 9.665, 25.98 and 

34.82 fim with the Short Wavelength Spectrometer. Atomic fine-structure lines were 
, detected from (in order of atomic number): C + , N + , N ++ , 0°, ++ , +++ , Si + , P + , 

and Fe + . The two lines of H2 that we observed, S(3) and S(9) were detected for all 
three remnants. The observations require both shocks into gas with moderate (~ 10 2 
cnr 3 ) and high (~ 10 4 cm" 3 ) pre-shock densities, with the moderate density shocks 
• producing the ionic lines and the high density shock producing the molecular lines. 

No single shock model can account for all of the observed lines, even at the order of 
magnitude level. We find that the principal coolants of radiative supernova shocks in 
moderate-density gas are the far-infrared continuum from dust grains surviving the 
shock, followed by collisionally-excited [O I] 63.2 /jm and [Si II] 34.8 /im lines. The 
principal coolant of the high-density shocks is collisionally-excited H2 rotational and 
ro-vibrational line emission. We systematically examine the ground-state fine structure 
of all cosmically abundant elements, to explain the presence or lack of all atomic fine 
lines in our spectra in terms of the atomic structure, interstellar abundances, and 
a moderate-density, partially-ionized plasma. The [P II] line at 60.6 //m is the first 
known astronomical detection, but its brightness can be explained using the solar 
abundance of P. There is only one, bright unidentified line in our spectra, at 74.26 
/um; as there is no plausible atomic fine-structure line at this wavelength, we suggest 
this line is molecular. The presence of bright [Si II] and [Fe II] lines requires partial 
destruction of the dust. The required gas-phase abundance of Fe suggests 15-30% of 
the Fe-bearing grains were destroyed. Adding the Si and Fe gas mass, and correcting 
for the mass of other elements normally found in dust, we find ~ O.5M of dust 
vapors from the shocked clump 3C 391:BML. The infrared continuum brightness 
requires ~ 1M Q of dust survives the shock, suggesting about 1/3 of the dust mass was 
destroyed, in agreement with the depletion estimate and with theoretical models for 
dust destruction. 

Subject headings: supernova remnants, shock waves, infrared: ISM: lines and 
bands, line: identification, ISM: abundances 
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1. Introduction 

Infrared spectroscopy is a particularly useful tool for studying supernova remnants (and 
other shocks) in dense gas for the following reasons. (1) Infrared atomic fine-structure lines are 
produced by gas with the physical conditions expected behind radiative shocks. Dense, shocked 
gas rapidly cools to temperatures of 10 2 -10 3 K, where infrared transitions are excited. This 
cool layer dominates the shocked column density, while the layer that emits optical to X-ray 
lines is extremely narrow. The kinetic temperature of the cooling post-shock gas is perfectly 
suited to the energy levels of infrared fine-structure lines from atoms and ions and infrared 
vibrational, ro- vibrational, and quadrupole vibrational (H2) lines from molecules. (2) Infrared 
fine-structure lines trace the ground-state populations of the all of the abundant elements that 
have electronic ground-states with non-zero angular momentum. By measuring the brightness of 
transitions directly to the ground-state or within a narrow multiplet of energy levels including 
the ground-state, a reasonably accurate measure of the column density of each species can be 
determined. Using pairs of lines from different elements due to energy levels with similar spacings, 
it is possible to measure relative abundances with confidence. (3) Far-infrared fine-structure lines 
are not affected by extinction. The light produced by post-shock gas is significantly affected 
by extinction from dust in the pre-shock cloud and the intervening medium along the line of 
sight. Extinction effects are severe in the ultraviolet and also in the optical and near-infrared for 
inner-galaxy supernova remnants, where the line-of-sight column density often exceeds 10 22 cm -2 , 
but it is negligible in the far-infrared. Many supernova remnants interacting with clouds are known 
only from their radio or hard X-ray emission, despite being copious producers of soft X-rays and 
optical line emission, because they are too distant or they are within or behind a molecular cloud. 

Recent advances in infrared astronomy have opened the far-infrared window to spectroscopy 
of faint emission. The observations used in this paper were made with the Long- Wavelength 
Spectrometer (LWS; plegg et al. 199(f ) aboard the Infrared Space Observatory (ISO; [Kessler et 



al. 1996|) , whose mission lasted from 1996-1997. A complete spectrum with the LWS, covering 



the 40-190 //m wavelength region with a resolution of 0.29-0.6 /im, took less than half an hour to 
reach a sensitivity adequate to detect atomic fine-structure lines and molecular rotational lines 
for a variety of astronomical sources. Instruments planned for the Stratospheric Observatory 



for Infrared Astronomy (SOFIA; Becklin 1997 ; expected to begin observing in 2002) will open 
opportunities for far-infrared spectroscopic observations over most of the far-infrared wavelength 
range — specifically, the part not affected by absorption and emission by very abundant molecules 
in the residual atmosphere such as H2O. In the somewhat-more-distant future, instruments aboard 
the Far-Infrared Space Telescope (FIRST; |Poglitsch 199£ ; anticipated launch 2007) will allow 



continuous spectral coverage to even deeper levels for extended sources. 

The goal of this paper is to present our ISO LWS and SWS spectral observations of the 
molecular supernova remnants W 44, W 28, and 3C 391, and to illustrate from basic physical 
principles why the spectral lines we detected are bright. Shocked gas in supernova remnants emits 
virtually all fine-structure lines from the abundant elements. We show the inferred elemental 
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abundances after impact from the supernova shock, and we discuss implication of grain destruction 
for elements such as Fe, Si, and P which were locked in grains. 

This paper continues our investigation of infrared emission from molecular supernova 
remnants. The brightness of the [O I] 63.2 /im line for W 44 and 3C 391 was presented in our 
first paper flReach fe Rho 1996| ), The far-infrared molecular emission of OH, CO, and H2O was 
presented in the second paper flReach fe Rho 1998 ). The millimeter- wave molecular emission of 



CO, CS, HCO + was presented in the third paper for 3C 391 QReach fe Rho 1999|) , and in this 



paper we also present some molecular results for W 44 from new IRAM 30-m observations. Two 
recent and closely-related studies of molecular emission were presented by Seta et al. (1998) for 
W 44 and Arikawa et al. (1999) for W 28. The millimeter-wave observations show beyond doubt 
that all three of the supernova remnants discussed in this paper are interacting with molecular 
clouds. 
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Fig. 1. — ISO LWS spectra of molecular shocks in the supernova remnants W 28, 3C 391, and 
W 44- The lower and upper panels show the wavelength ranges 42-80 and 80-187 \xm, respectively. 
The ordinate in both panels is vl v ; to measure line fluxes note that the resolution of the LWS is 
= 0.29 \xm for the top panel and 0.6 fim for the bottom panel. Note that no reference spectrum 
has been subtracted, so the continuum and [CllJ line are partly due to line-of-sight dust and gas. 
The other bright atomic fine structure lines are due to the remnants. Some of the minor wiggles at 
long wavelengths are molecular lines (Reach & Rho 1998), while detector noise is significant at the 
shortest wavelengths. The unidentified line is labeled with a '?'. 
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2. Observations 

2.1. Far-infrared spectral lines 

We obtained complete ISO LWS spectra of 1 position in each of three supernova remnants as part 
of our program of observing supernova-molecular cloud interactions (Reach & Rho 1996, 1998, 
1999). The positions observed were the brightest OH masers ( Frail et al. 1996|) in W 28 



(18 fc 01 m 52.3 s -23°19'25") and W 44 (18 h 56 m 28.4 s +1°29'59") and the molecular peak in 
3C 391:BML (18 /l 49 m 21.9 s -0°57'22"); all coordinates are J2000. Figure g shows the ISO LWS 
spectrum from 42-188 fan for all three remnants. For each of the 10 LWS detectors, we removed 
the fringes using the ISO Spectral Analysis Package^] and subtracted a constant (due to residual 
dark current). Table [l] shows a list of the identified atomic fine structure lines and their measured 
brightness. In order of atomic number, we detected atomic fine-structure lines from C + , N + , 
N++, 0°, 0++, 0+++, and P+ in the LWS spectra. 

Part of the bright continuum and [C II] 157.7 /jm line in Fig. [l] are due to dust and gas along the 
line of sight, unrelated to the supernova remnant. At least part of the continuum is likely to be 



due to dust grains surviving the shocks. For example, for W 44 we estimated 30% ( Reach fc Rho 



1996) based on comparison to a reference position; a spectral analysis of the continuum is given in 
§H of this paper. To estimate the line-of-sight contamination of the spectral line brightnesses, we 
calculated the galactic line-to-continuum ratio using the COBE Far-Infrared Absolute 
Spectrophotometer spectrum of the galactic plane at longitude 30° ( [Reach et al. 1995| ). At the 
ISO LWS resolution, the expected line-to-continuum ratio due to unrelated gas for [C II] 157.7 //m 
is 2.2, while that for [N II] at 121.9 fim is 0.05. Comparing to our spectra, the 157.7 /im line to 
continuum ratios are 2.4, 1.2, and 1.9, and the 121.9 /Lim line to continuum ratios are 0.19, 0.09, 
and 0.16, for W 44, W 28, and 3C 391, respectively. Thus much of the [C II] /jm line emission 
could be due to line-of-sight gas, while only a small fraction the [N II] (and essentially none of the 
other lines) are due to line-of-sight gas. The subject of [C II] emission being related to the shock 
is addressed again in §|4|, focusing on the abundance of C. An improved separation of the source 
and background brightness might be possible by analyzing the spatial variation of the brightness, 
which is beyond the scope of this paper. 



The ISO Spectral Analysis Package (ISAP) is a joint development by the LWS and SWS Instrument Teams and 
Data Centers. Contributing institutes are CESR, IAS, IPAC, MPE, RAL and SRON. 
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Fig. 2. — ISO SWS spectra of the ground-state [Fe 11/ fine-structure line for lines of sight toward 
molecular shocks in (a) W 28, (b) W 44> an d ( c ) 391; also the spectrum for a line of sight 
toward the radio peak in 3C 391. The molecular shock positions all show bright [Fe II 7 lines. The 
radio peak position in 3C 391 also shows a bright ground-state fune structure line of [O TV]. The 
solid curves are Gaussian fits to the line with width fixed to the SWS grating spectral resolution of 
AX = 0.018 fim. 



- 7- 




Fig. 3. — ISO SWS spectra of the ground-state [Si II/ fine-structure line for lines of sight toward 
molecular shocks in (a) W 28, (b) W 44> an d (c) 3C 391; also the spectrum for a line of sight 
toward the radio peak in 3C 391. The solid curves are Gaussian fits to the line with width fixed to 
the SWS grating spectral resolution of AA = 0.022 fim. 
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Fig. 4. — ISO SWS spectra of the H2 S(9) rotational line toward molecular shocks in (a) W 28, (b) 
W 44: an d ( c ) 3C 391; also the spectrum for a line of sight toward the radio peak in 3C 391. The 
solid curves are Gaussian fits to the line with width fixed to the SWS grating spectral resolution of 
AA = 0.0021 [im. 
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Fig. 5. — ISO SWS spectra of the H2 S(3) rotational line toward molecular shocks in (a) W 28, (b) 
W 44: an d ( c ) 3C 391; also the spectrum for a line of sight toward the radio peak in 3C 391. The 
solid curves are Gaussian fits to the line with width fixed to the SWS grating spectral resolution of 
AA = 0.0042 [im. 
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Table 1: Observed far-infrared Atomic fine-structure lines 



Line Brightness (10 erg s cm sr ) 



ion transition wavelength 



C+ 

0° 

N+ 

0++ 

0° 

p+ 

N++ 
0++ 



W 28 



W 44 



3c 391:BML 



2p 



3 Po^i 
3 P 2 -^i 
3 Pi^o 

3 Pl-2 
3 Pl^0 

3p 



l 

'T2 



2->l 



157.741 
145.525 
121.898 
88.356 
63.184 
60.640 
57.340 
51.815 



7.77 ±0.08 
0.75 ±0.04 
0.68 ± 0.04 
0.38 ±0.06 
9.48 ± 0.02 
0.06 ±0.03 
0.093 ± 0.033 
0.23 ±0.06 



4.94 ± 0.04 

0.50 ±0.03 
0.42 ± 0.04 
0.36 ± 0.02 
7.73 ± 0.04 

< .10 

< .14 
< 0.23 



6.54 ±0.08 
0.90 ±0.02 
0.89 ± 0.02 
1.49 ±0.04 
12.96 ±0.02 
0.095 ± 0.025 
0.25 ±0.04 

0.75 ±0.08 



Table 2: SWS Spectral line results 



ion 


transition 


wavelength 


Line Brightness (10 


4 erg s 1 cvci 


^sr- 1 ) 


W 28 


W 44 


3c 39LBML 


3C 391:radio 


Fe+ 


2 2 


25.988 


0.66 ± 0.05 


0.41 ± 0.04 


1.9 ±0.1 


1.8 ±0.1 


o+++ 


2 P3^! 

2 2 


25.890 


< 0.08 


< 0.08 


< 0.08 


0.12: 


Si+ 




34.815 


7.3 ±0.1 


4.8 ±0.1 


15 ±0.2 


13 ±0.2 


H 2 


2 2 

S(3) 


9.665 


10.2 ±0.1 


3.6 ±0.1 


1.95 ±0.05 


0.32 ± 0.06 


H 2 


S(9) 


4.695 


1.2 ±0.1 


1.9 ±0.1 


1.3 ±0.1 


< 0.2 
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2.2. Mid-infrared Si, Fe, and H2 lines 



We used the ISO SWS (de Graauw et al 1996) to observe narrow ranges centered 25.98 and 34.82 
/jm for the ground-state transitions of [Fe II] and [Si II]. The spectra are shown in Figures ^ 
and HI respectively. At the same time, we obtained spectra of the H2 S(9) and S(3) lines; the 
spectra of these lines are shown in Figures || and || respectively. The brightnesses of the lines are 
summarized in Table [2|. The SWS observations were toward the three lines of sight listed above, 
plus one extra position — toward the bright, northwestern radio ridge and a bright peak in the 



[O I] 63 lira strip map (position b in Fig. 3 of [Reach k Rho 1996| ): 3C 391:radio (18 /l 49 m 19.4 s 



-0°55'05"). The two positions in 3C 391 have similar spectra of ionic fine structure lines: the ratio 
of [Fe II] /[Si II] is nearly identical, although the [O IV] is relatively brighter toward 3C 391:radio 
by at least a factor of 3. In contrast, the brightness of the H2 lines, relative to the ionic lines, is 
more than a factor of 10 lower for 3C 391:radio than for 3C 391:BML. This is consistent with the 
millimeter-wave molecular lines revealing no evidence for shocked molecular gas toward 
3C 391:radio, while bright and broad shocked molecular lines were detected toward 3C 391:BML 
(Figs. 2 and 4 of Reach fc Rho 1999| ). The spectra of these two positions serve a valuable 



comparison and diagnostic tool, because one of them is dominated by ionic lines, while the other 
has significant molecular emission. 
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Fig. 6.— (a) High-resolution ISO LWS spectrum of the O I line toward 3C 391:BML. The 
instrumental resolution is 36 km s^ 1 , so the line is clearly resolved, (b) and (c) Spectra of the 
CO (2 — > I) and 13 C0(1 — > 0) lines for the same line of sight, obtained with the IRAM 30-m 
telescope. The 13 CO line arises only from unshocked gas, while the CO (2 — ► 1) line has a wide 
component from the shocked gas as well as narrow components from the unshocked gas. 
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Fig. 7. — (a) High-resolution ISO LWS spectrum of the Ol line toward W44- (b) and (c) Spectra 
of the CO(2 — > 1) and 13 CO(l — > 0) line for the same line of sight, obtained with the IRAM 30-m 
telescope. 
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2.3. High-resolution [O I] and CO spectra 

We used the Fabry-Perot in the LWS to obtain high-resolution spectral profiles of the [O I] 63 fim 
lines for 3C 391:BML and W 44. The positions are the same as those where the complete LWS 
spectra and the SWS line spectra were taken, and they correspond to peaks in the [O I] strip 
maps presented in paper 10. The instrumental resolution of the LWS Fabry-Perot is 36 km s _1 . 
The two spectra are shown in Figures ^ and fjj. The lines are clearly resolved, with a full width at 
half maximum (FWHM) of 100 km s" 1 . The lower panels of Figs. | and @ show the CO(2 -> 1) 
and 13 CO(l — > 0) line profiles for the same lines of sight as the O I spectra, obtained with the 
IRAM 30-m telescope ( [Wild 1995 ). The observing techniques are described in our previous paper, 



where the 3C 391 observations were presented (Reach & Rho 199£); the W 44 observations are 



new. The line profiles are all clearly, very different. The 13 CO lines have a narrow (1.5 km s _1 
FWHM) velocity dispersion, and they trace the column density, so they arises from large clouds 
of cold, unshocked or unrelated gas. Although the CO(2 — > 1) lines also have narrow components, 
identical to the 13 CO lines, from the unshocked or unrelated gas, the CO(2 — > 1) lines are 
dominated by a broad component (20-50 km s _1 FWHM). The broad component was seen in the 
CO, CS, and HCO + observations of 3C 391, with relatively more emission from the 



higher-excitation and higher dipole moment lines ( Reach fc Rho 1999 ). Similarly, the CO 
observations of W 28 showed that the broad component was much brighter than the narrow 
component in the 3 — > 2 line as compared to the 1 — > line ( Arikawa et al. 1999| ). The broad 



components of the molecular lines are still not as broad as the [O I] line, suggesting that each of 
the spectral lines in Figs. ^ and [7| trace different regions. 

Inspecting the [O I] line profiles, there is only slight evidence for structure. The lines from both 
remnants can be reasonably represented by gaussians with FWHM ~ 100 km s _1 . For 
3C 39LBML, the line profile is somewhat better represented by multiple components: an 
unresolved component centered near 50 km s _1 and a broad, non-Gaussian component centered 
near km s _1 with 120-150 km s -1 . The integrated line brightness is dominated by the wider 
component. The unresolved component contributes only about 10% of the line integral, and it 
may be associated with preshock or unrelated gas, which would agree with the brightness of the 



'off-remnant' positions from the [O I] strip map ( Reach fc Rho 1996 ). Neither 3C 391 nor W 44 



show evidence for self-absorption in the line profiles. This was a particular concern in interpreting 
the line brightnesses for molecular supernova remnants, both because of the possibility of dense 
pre-shock clouds and the large column densities of cold gas along the lines of sight to the 
remnants. A single cold cloud might be expected to yield an essentially unresolved (Gaussian) 
absorption dip, which does not occur. The foreground gas would be expected to absorb in 



2 Please note that the apparent 400 km s _1 shift in line center in 3C 391 from one side of the remnant to the 
other, mentioned in paper 1, is apparently an instrumental artifact due to mapping in the dispersion direction of the 
grating. The observations presented in this paper, with the Fabry-Perot, are considered reliable measurements of the 
line profile. 
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approximate proportion to the line of sight column density. After inspecting the H I 21-cm profile 
for these lines of sight, we found no correlation that could be due to absorption dips. Instead, it is 
likely that the line-of-sight gas contributes slightly to the emission. The level of [O I] emission 
from the line-of-sight gas is most easily estimated from the reference positions just outside the 
remnant, where it was found to be quite faint ( [Reach fe Rho 1996j ). Therefore, we interpret the 
[O I] line profiles in Figs. [6| and ^ as emission from the shocked gas. This inference yields a direct 
measure of the shock velocity of 100-150 km s _1 , albeit in projection along the line of sight. 
There are likely to be multiple shocks within the beam, and the FWHM of the line profile is 
approximately equal to the 'typical' shock velocity. 

Comparing the [O I] and CO line profiles, we can partially describe the three-dimensional 
geometry of the shocks. For 3C 391, note that the centroids of the [O I] and CO lines are very 
different. The CO line occurs at nearly the systemic velocity at the remnant's distance, with the 
broad component due to the shocked gas very slightly blueshifted from the narrow line from the 
preshock gas. The [O I] line, on the other hand, is significantly red-shifted from the pre-shock gas. 
Thus the shock that produces the [O I] line has a significant velocity component toward us; that 
is, the shock is propagating partially along the line of sight. This suggests that the 3C 391 
progenitor exploded on the far side of its parent molecular cloud, driving the dense blast wave 
partially in our direction. For W 44, the centroids of the [O I] and CO lines are similar, 
suggesting that the shocks are either perpendicular to the line of sight, or that the lines come 
from similar geometries. One straightforward conclusion from comparing the line profiles, both in 
centroid and width, is that the [O I] line does not arise exclusively from slow shocks into very 



dense gas such as produces the CS, OH, and H2O emission ( Reach Rho 1999 ). Instead, most of 



the [O I] emission must arise from faster shocks into moderate-density gas. 



3. Comparison of detected lines to the periodic table 

3.1. Atomic fine-structure lines: Basic principles 

The three basic principles that determine which spectral lines will be bright in the spectrum of a 
parcel of gas are abundance, ionization, and excitation. The abundances of the elements are given 
by Anders & Grevesse (1989) for the Solar System. In the interstellar medium, since part of the 
abundance will be locked in solids, the gas phase abundance will be smaller than Solar System 
abundance for some elements ( [Savage fc Scmbach 1996] ) . It is clear that we need only consider the 
top three rows of the Periodic Table and the middle of the fourth row; all other elements have 
negligible abundance. The ionization state of each element depends on where the gas is located. 
We can consider three typical locations: the diffuse interstellar medium (where all electrons bound 
by less than 13.6 eV are removed), dark clouds (where elements are neutral), and hot plasmas 
(where the higher ionization states are present). The excitation of the atomic levels will depend 
on the density and temperature of the region, but a 'zero-th' order effect determines whether 
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there are any atomic fine structure lines in the ground electronic state. Elements and ions with 
zero angular momentum {L = 0) in their ground state electron configuration will not have the 
fine-structure lines considered here. For example, in spectroscopic notation, S configurations have 
no lines, while doublet 2 P configurations have 1 line and triplet 3 P configurations have 2 lines. 
The fine-structure lines are magnetic dipole transitions, so they are technically 'forbidden lines', 
but with spontaneous radiative decay rates of order 10~ 5 s _1 and collisional de-excitation rates of 
order 10~ 7 cm -3 s^ 1 , they are easily excited under astronomical conditions and are thermalized 
at densities ;> 10 2 cm -3 (for far-infrared lines) to 10 4 cm -3 (some mid-infrared lines). 

In order to understand and predict which fine-structure lines are present in astronomical spectra, 
we constructed a 'periodic table for astronomical fine structure lines'. Table || contains the 
electron configurations and wavelengths for the p-shell ions of the second and third rows of the 
Periodic Table (C, N, O, Fe, Ne, Al, Si, P, S, CI, and Ar), while Table || has the same information 
for the fourth row of the Periodic Table (Cr, Mn, Fe, CO, and Ni). In both tables, the dominant 
ionization state in the diffuse interstellar medium is shown in bold. Spectra of regions behind 
shock fronts or with very high radiation fields will contain the dominant ionization state as well as 
and the ions to the left of the dominant ionization state. Spectra of regions that are shielded from 
the interstellar radiation field will contain neutral forms and ions to the right of the dominant 
ionization state. Some elements with negligible astronomical abundance are not listed; also, to 
reduce clutter, unlikely ionization states of rare elements are not listed. Alkaline metals (in low 
ionization states) have no p-shell electrons to make infrared fine-structure lines, so they are not 
listed. These tables can be used both to predict which lines will be present for a given gas and to 
decide which lines should be present when another line of the same ion is detected. In all cases, 
the transitions are simple cascades down a ladder, with the ground state at the bottom; therefore, 
the upper energy level for a line can be calculated by summing the photon energies from a given 
line to the bottom. 



3.2. Comparison of detected lines to the Periodic Table 

In the 'periodic table for atomic fine-structure lines' (Tables [| and Q), the lines we detected from 
molecular supernova remnants observations are in italics. It appears that we detect all of the 
available ground-state fine-structure lines of the astronomically abundant molecules; the 
undetected ground-state lines are generally outside our wavelength range. All abundant elements, 
C, N, O and Si, are detected as shown in Tables |l] and |2[ Non-detections of F and Co lines are 
due to low abundances. Lines from high- ionization states were also not detected. The only 
element with high abundance that we did not detect is Ne, because our observations did not cover 
the wavelengths of its ground-state fine-structure lines. Our own observations contain the 
wavelengths of 12 fine-structure transitions among energy levels within 600 K of the ground state, 
from elements with abundance > 10 -7 relative to H in their expected ionization states. We detect 
essentially all of these lines, as well as lines from higher ionization states of the abundant 
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Table 3. 'Periodic Table' for Astronomical Fine Structure Lines: Z=6— 18 



Element 


Abundance 13 


2 Pl/2 


3 Po 


4 S 3 /2 


3 P 2 


2 P S /2 


1 Ss/2 




C 


-3.4 


Ul57.7 


,370 
l 610 










N 


-3.9 


Ulsr.s 


rrlB1.9 
ii 205.9 


I... 








O 


-3.1 


IV 25.9 


TTT5J.S 


II... 


X 6S.2 






F 


-7.5 








T 67.2 
i 29.3 


I24.8 




Ne 


-3.9 






IV ... 


rrr36.0 
JII 15.6 


Hl2.8 


I... 


Al 


-5.5 


189.2 












Si 


-4.4 


U34.8 


T 68.4 
J-129.7 










P 


-6.4 


III17.9 


TT 32.9 
11 60.0 


I... 








s 


-4.7 


IVlO.5 


TTT18.7 

m 33.5 


II... 


T 56.3 
J-25.2 






CI 


-6.7 






III... 


TT33.3 
li 14.4 


111. 3 




Ar 


-5.4 






IV... 


TTT21.8 


He. 99 


I ... 



a The headings of columns 3-8 give the spectroscopic notation X Y Z for the ground 
state. Roman numerals represent the ionization state (I=neutral, II=singly ionized, 
...), with the dominant ionization state in the diffuse interstellar medium in bold. 
Small numbers indicate the wavelength of the ground-state fine structure lines, in 
fim. The transition to the true ground state is at the bottom and increasing energy 
levels are above, when applicable. All wavelengths of spectral lines that we detected 
are shown in italics. Zero-spin ground states, which have no fine structure lines, have 
in their wavelength column. 

h A is the base-10 logarithm of the abundance relative to H. 



Table 4. 'Periodic Table' for astronomical Fine Structure Lines: Z=24-28 



Element 


Abundance 


23 


24 


25 




26 


27 


28 


Cr 


6.3 


II 6 S 5/2 ... 


i 7 s 3 ... 












Mn 


6.5 


IH 6 S 5 /2... 


II 7 S 3 ... 


I 6 S5/2 










Fe 


4.5 


IV 6 S 5/2 ... 


1U5.4 

III 5 D 4 ^; 7 
22.9 


II 6 D 9/2 


87.3 
51. 3 
35.3 
25.99 


111.2 

t 5r-v 54.3 
24.0 






Co 


7.1 






III 4 F 9/2 


24.1 
16.4 
11.9 


Tt3t? 15.5 
11 i4 10.5 


t4tt< 16.9 
1 * 9/212.3 




Ni 


5.7 










IH 3 F4 £S 


H 2 D5/26.64 


T 3-p 11.3 
1 47.51 



a Column heading is the number of electrons. See caption for Tab. 
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elements: O III, O IV, and N III. In addition to these, the ground-state line of [P II] is also detected 
for W28 and 3C391. Some of the other lines were detected in the younger and brighter supernova 
remnant RCW 103 flOliva et al. 1999bfl ; upon reanalysis of these data from the ISO archive, we 
find that the detected lines include the ground-state line of [P II]. We can use these general 
principles to predict other bright lines, outside of our observed spectral range. For example, Ne is 
a very abundant element with lines in the mid-infrared; given that we detected ionized states of 
O, we expect bright emission from [Ne II] at 12.8 fim and perhaps [Ne III] at 15.5 /um. 



10 6 ^ 



10 5 - 



^ m 4 " 
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10 3 - 
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region too large 




100 1000 

Temperature (K) 



10000 



Fig. 8. — Excitation diagrams for the [Ol] (upper panel) and [Olll] (lower panel) lines. The solid 
contours show the observed range of 145/63 \im (upper panel) and 52/88 jj,m (lower panel) line 
ratios from Tabel [/j. 
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The dashed contour running approximately from the upper left to lower right of each panel shows 
the constraint placed on the density and temperature by requiring that the emitting region is 
smaller than the entire supernova remnant, while still producing lines as bright as observed. The 
vertical dashed line in the upper panel at T = 10 4 K shows where O should be collisionally ionized. 
The dashed contour near the top of the lower panel shows the constraint that the O III region be at 
least as hot and no more dense than the O I region. Hatched regions in both panels show the 
allowed solutions. The asterisk in the upper panel shows the solution that we chose for detailed 
calculations. The upper-left asterisk in the lower panel shows this same solution if the gas is fully 
ionized, and the asterisks connect to it show the set of solutions in thermal pressure equilibrium 
with it. The lower set of asterisks shows the same solutions assuming the ionization fraction in 
the O III region is 1%. 



4. Abundances of the shocked gas 

4.1. Excitation of Oxygen 

Based on the arrangement of the energy levels of the ground-state atomic fine structure, the 
ratios the upper and lower [O I] and [O III] lines constrain the density and temperature of the 
emitting region, independent of the geometry. The upper panel of Figure || shows the line ratio 
predicted for [O I] as a function of the total H density and temperature. At each density and 
temperature, the level populations were calculated and the emergent intensity predicted using the 
escape probability ( Hollenbach fc McKee 1979| ). The column density, Ah, was adjusted iteratively 



for each model such that the model predicts the correct brightness of the [O I] 63 /im line. Solid 
lines show the locus of models that predict the observed line ratios from Table |]. The locus of 
models for which the path length through the emitting region, z = Nn/nn, is equal to 3 pc (i.e. 
comparable to the radius of the supernova remnant) is shown as a dashed contour, and the region 
below this contour is labeled 'region to large.' Low-temperature and low-density models are not 
capable of producing the observed line brightness within the confines of the observed region. 
Another constraint is drawn as a vertical line at T = 10 4 K, where the neutral O atom would be 
collisionally ionized; higher temperature models should therefore be excluded. The lower panel of 
Figure |8| shows the line ratio predicted for [O III] as a function of density and temperature, with 
each model normalized to predict the correct brightness of the [O III] 88 /im line. The [O III] and 
[O I] emitting regions are not necessarily coincident, so we consider the density, temperature, and 
ionization in the two regions separately. One relative constraint between the two emitting regions 
is practically guaranteed: the temperature will be at least as high, the density will be at least as 
low (and the ionization will be at least as high) in the [O III] emitting region as compared to the 
[O I] emitting region. This rules out the very dense models for the [O III] emitting region, as 
indicated by the upper dashed curve in the lower panel of Figure 0. 



To determine the physical conditions where the lines are produced, and to determine the 
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abundances of the elements in the shocked gas, we define two emitting regions. Region M 
(molecular) is the [O I] emitting region, with density, temperature, and column density set to give 
the observed brightnesses of the 63 and 145 /im lines. Region A (atomic) is the [O III] emitting 
region, with density, temperature, and column density set to give the observed brightnesses of the 
52 and 88 /im lines. Such a separation of emitting regions is anticipated because we expect that 
the neutral and doubly-ionized states of O would not be co-spatial — although there could in 
principle be overlap, if the gas is far from ionization equilibrium. Combining all of the 
constraints, the plausible range of density and temperature for region M has density 
2000 < n^ 7) (cm~ 3 )< 5000 and temperature 100 < (K)< 1000 K. For convenience, we have 

chosen a nominal solution of Wjj = 3000 cm -3 and T^ M ^ = 600 K, which is marked with an 
asterisk (*) in Figure || and listed in Table ||, for detailed calculations. 

In region A, there is a wider range of solutions within the constraints already mentioned. In 
particular, models with constant n e = n^x, where x is the fraction of H that is ionized, produce 
the same excitation for the ions in region A. If region A and region M were layers behind the 
same shock front, then we expect region A to have a comparable or somewhat higher pressure 
than region M (lest the shock move backwards). For the nominal solution for region M, the region 
A solutions with the same pressure but density lower by factors of 1, 2, 4, and 8 are shown as two 
lines connecting asterisks in the lower panel of Figure ||. The upper line applies if region A is 
completely ionized, and the upper left asterisk is just a direct copy of the region M solution. It is 
evident that region A cannot be completely ionized and be a shock layer upstream from region M: 
the region A pressure would be so low that the shock would move backwards. The lower line 
applies if region A has an ionization fraction x = 10 -2 , and the upper- left asterisk of the lower 
line is the region M solution with n e = nn /W0. This would be a very surprising result, because it 
means that ++ coexists with a significant amount of H°. Very low density gas can be far from 
ionization equilibrium, with higher ionization states coexisting with neutral states (Schmutzler & 
Tscharnuter 1993), but at the high densities required to get the observed brightness of the 63 and 
88 fim lines the ionization should be reasonably close to equilibrium. For this reason, we consider 
the [O III] emitting region (region A) separate from the [O I] emitting region. For detailed 
calculations, the region A model has density = 8 cm -3 and temperature = 3000 K; the 
properties are listed in Table || In region A, the abundance of O is set to zero, and the 
abundance of ++ and N ++ are set to the entire cosmic abundance of O and N, respectively. 



4.2. Predicted line brightnesses for different regions 

The impact of a supernova blast wave with a molecular cloud leads to a complex interaction, with 
regions of different density reacting in very different ways. Based on the excitation of [O I] and 
[O III] (discussed in the previous section), and the excitation of H2O and OH ( [Reach fc Rho| 
|1998[ ), we identify three types of post-shock gas and their infrared emitting lines. Figure |9]is a 
cartoon showing these shocks, and Table [| lists their properties. The lowest-density pre-shock gas 
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is region A (for 'atomic') with pre-shock density Uq < 1 cm -3 . (If the pre-shock density of 
region A were higher than 1 cm -3 , the post-shock density would be much higher than the 8 cm~ 3 
and would be inconsistent with the 52/88 (j,m. line ratio.) The post-shock gas in region A gives 
rise to bright H II recombination lines and infrared and optical O III and N III (and other high 
ion) forbidden lines. The moderate-density pre-shock gas is region M (for 'molecular') with 
pre-shock density tiq ~ 10 2 cm~ 3 . The post-shock gas in region M gives rise to bright O I and 
most of the other bright infrared forbidden lines we observed, as well as some optical emission 
(which is, however, highly obscured). The highest-density pre-shock gas is in the dense clumps, 
with pre-shock density Uq ~ 10 4 cm~ 3 . The post-shock gas in region C gives rise to bright H2, 
OH, CO, CS, and H2O emission, with a possible contribution to the O I emission. We do not 
mean to imply that there are three sharply-defined regions with fixed densities; we do suggest 
that the full range of pre-shock densities is needed, such that no single region is consistent with 
all of the observed spectral lines. 

The shock velocities in the different regions are not well known, but reasonable estimates can be 
made. If region A is the lowest-density pre-shock gas, then it should represent the leading edge of 
the supernova remnant. The Sedov solutions for these remnants are consistent with the density of 
region A for explosions with energy of order 10 51 erg ( [Rho 1995 ). Thus the shock velocity should 



be of order 2R/5t, where R is the radius and t is the age of the remnant. This works out to 
vj A ^ ~ 500-600 km s" 1 for the remnants considered in this paper ( [Rho 1995|) . For region M, we 



have some direct evidence of the shock velocity from the width of the O I line (§ |2.3|) that 
y( M ) iqq k m g -i_ This i s a i so consistent with the expanding H I shell observed toward W 44 



( Koo &: Heiles 1995| ). For region C, there is some direct evidence of the shock velocity from the 
width of the millimeter-wave CO and CS lines flReach fc Rho 199GQ that V s { ' ~ 30 km s" 1 . The 



shock ram pressure p ram ot noV^ should be comparable in these regions as they are all driven by 
the same blast wave. Indeed we find p ram is progressively higher for the denser regions, but the 
range of p ram is only a factor of 36. A pressure enhancement in the denser regions has been 
theoretically justified by Chevalier (1999) as due to the shock into the denser gas (regions M and 
C) being driven from the radiative shell (region A) as opposed to being isolated shocks driven by 
the blast wave. 
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RegionA, ::: .>.„. Region M, 




H 2 , H 2 0, OH 



Fig. 9. — Cartoon illustrating the different pre- and post-shock regions discussed in the paper. The 
physical properties of the regions are listed in Table |^. For each region, we also list an ion or 
molecule whose emission can be used to image the region. 

In Chevalier's (1999) paper, two regimes of density were considered for the pre-shock gas; these 
densities correspond roughly to our regions A and M (cf . Table [BJ) . Chevalier explained the [O I] 
63 /im luminosity of the remnant using only relatively low-density shocks into region A. Based on 
our new results for the line ratio and excitation of O I and the high-resolution spectrum of the 63 
//m line showing a width of only 100 km s _1 , we suspect that the 63 //m arises in relatively denser 
regions, like region M. Indeed, it is impossible for a low-density emitting region to produce 63 /xm 
lines as bright as we observed, within the confines of the supernova remnant. In the model of 



Table 5: Physical properties of the three post-shock emitting regions 





pre- 


shock 


post-shock 






Region 


n 


V s 


n 


T 




N 


Infrared Coolants 




(cm -3 ) 


(km s _1 ) 


(cm -3 ) 


(K) 


(' 


cm -2 ) 




A (atomic) 


< 1 


500 


8 


3000 




10 20 


O III, N III 


M (molecular) 


10 2 


100 


3 x 10 3 


600 




10 21 


O I, Si II,Fe II, N II, H 2 S(9) 


C (clump) 


10 4 


20 


4 x 10 5 


200 




10 22 


CO, H 2 0, OH, H 2 S(3) 



-23- 



Hollenbach & McKee (1989), the 63 /im line brightness is proportional to uqV s . For the cartoon 
model parameters in Table ||, it is clear that the denser regions will dominate the line brightness if 
they fill an appreciable portion of the beam. In region C, however, the O is probably largely tied 
up in molecular form; furthermore, region C is somewhat beam diluted in the ISO LWS beam. 
Therefore, we believe that the 63 /im lines presented here arise from shocks into the 
moderate-density molecular gas of region M. 

For regions A and M, we calculated the model brightness of all of the lines that we observed with 
ISO. Each element was assumed to have its entire solar abundance in the gas phase and in the 
ionization state that gives rise to the line, with the following exceptions. In region M, oxygen was 
assumed to be neutral (as we observed), and nitrogen was assumed to be singly ionized (since the 
only detected line was from N II, not N I) in order to obtain an upper limit to the predicted N II 
brightness.. In region A, oxygen was assumed to be ionized. Our purpose in these calculations is 
to measure the richness or depletion relative to solar abundances. Our calculations do not yield 
absolute abundances, because in reality the total abundance of each element is actually divided 
among the ionization states that are present. In the next section, we go through the list of 
elements and discuss their expected ionization state and the inferred total abundances in the 
post-shock gas. Table |6] shows the brightness predictions for regions A and M. We can conclude 
the following, without considering the ionization balance. First, it is possible to reproduce all of 
the observed fine-structure lines, except N III and O III using region M alone. (In fact, it is 
possible to get the N III and O III lines, if we could accept the multiply ionized species existing in 
region M, where H is mostly neutral.) Therefore, we expect to be able to derive abundances in 
the 'region M' post-shock gas. For region C, we discuss details in §|| below. 



4.3. Post-shock gas abundances 

We now step through each element and derive the total abundance of that element in the 
post-shock gas. Table |7| summarizes the abundances inferred for the post-shock gas. For each 
element X, [X/H] & is the solar abundance ( |Anders &: Grevesse 1989| ) and [X/H] b s is the total 
abundance of that element from our observations and model. 

Carbon — The only carbon line we observed was from C II, which is likely to be the dominant 
ionization state in region M. The observed line brightness compared to the region M prediction 
(Table |6|) requires essentially all of the carbon to be in the gas phase. However, a significant 
fraction of the 157 [im line could come from line-of-sight gas, so this result is inconclusive. The 
carbon abundance in the post-shock gas could be better addressed if we knew the brightness of 
the line just outside of the remnant. Using our current knowledge, we can only say that the ISO 
observations are consistent with complete destruction of the C-bearing dust grains if all of the 
emission comes from the remnant, or partial or no destruction if less than half of the emission 
comes from the remnant. In Table |?] we list the inferred depletion of carbon in the post-shock gas 
assuming all of the 157 //m line comes from the remnant; this yields an upper limit to the 



Table 6: Model predictions for fine-structure lines 



Ion Abundance line wavelength Atomic Observed Observed/Predicted 

(/im) Data a (erg cm -2 s" 1 sr _1 ) layer M 6 layer A c 



CII 


-3.44 


157.7 


1 


5.8 x 10~ 4 


1.2 




Nil 


-3.95 


121.9 


1 


2.7 x 10~ 4 


0.3 


10 


NIII 


-3.95 


57.3 


2,3 


9.7 x 10~ 3 




0.3 


01 


-3.07 


63.2 


4,5 


1.3 x 10~ 3 


1.0 




01 


-3.07 


145.5 


4,5 


9.1 x 10~ 5 


1.0 




III 


-3.07 


88.4 


1 


5.2 x 10~ 3 




1.0 


III 


-3.07 


51.8 


1 


2.8 x 10~ 3 




0.9 


Sill 


-4.45 


34.8 


5 


3.0 x 10~ 3 


0.5 


10 


P II 


-6.55 


60.6 


6,7 


8.6 x 10~ 6 


1.1 




Fe II 


-4.33 


26.0 


8 


2.7 x 10" 4 


0.3 


10 



"References for atomic data: 1 ( Spitzer 1978 ), 2 ( Froese Fischer 1983j ), 3 ( [Blum fe Pradhan 1992 ), 
4 (|LaunayfeRoueff 19771 ), 5 QHollenbach fe McKee 1989p , 6 flMendoza fe Zeippen 19821) , 7 ( [Krueger 
fe Zsyzak 1970|) , 8 flNussbaumer fe Storey 1982]) 

fe layer M model parameters: n = 3000 cm -3 , T = 600 K, O III and N III abundances set to zero, 
column density normalized to produce observed [O I] lines 

c layer A model parameters: n e = 8 cm -3 , T = 3000 K, O I abundance set to zero, column density 
normalized to produce observed [O III] lines 



Table 7: Abundances 
Element log [X/H] Q [X / H] obs / [X / H] Q 



C -3.44 < 1.1 

N -3.95 > 0.3 a 

O -3.07 1 

Si -4.45 0.5 

P -6.55 1.1 

Fe -4.33 0.3 



abundance of N is a lower limit because the dominant ionization state (N I) was not observed 
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gas-phase abundance, and therefore lower limit to the depletion. 

Nitrogen — We observed nitrogen in its singly- and doubly-ionized states. The dominant state in 
region M should be neutral, because the ionization potential of N I is greater than 13.6 eV. 
Comparing the observed to predicted brightness (Table 0), we require only about 1/3 of the 
nitrogen to be N II, which is reasonable. When all nitrogen is assumed to be singly ionized, the 
observed to predicted brightness ratio of [N II] is only 0.3 (Table ||), which implies either nitrogen 
has lower abundance or the predicted N II is too high. Since the former is unlikely, the ratio 
implies that only about 1/3 of the nitrogen is N II and the rest is neutral. In region A, which was 
normalized to give the correct 88 /jm line brightness, we require about 1/3 of the N III, which is 
also plausible. In Table we list the depletion of nitrogen assuming all of the nitrogen is ionized; 
this is a lower limit to the gas-phase abundance, and an upper limit to the depletion, because it 
neglects the neutral state. 

Oxygen — The column density of our model was adjusted to match the observed brightness of the 
oxygen lines, so most of our measurements are, effectively, relative to the O abundance. We 
cannot easily determine the total abundance of O because we lack an accurate measurement of 
the total column density of H. 

Silicon and Iron — We find that a significant fraction of the total abundance of Si and Fe, which 
are normally highly depleted into dust grains in the interstellar medium, is required to produce 
the observed line brightness. Both Si and Fe were observed in their dominant ionization state in 
region M. Region A does not have enough column density to contribute significantly to the 
observed line brightness. (Recall that the region A column density is set by the [O III] line 
brightness, so increasing its column density would make too much [O III] emission.) For region M, 
we find that about 1/2 of the Si and 1/4 of the Fe are required to be in the gas phase. These 
abundances of Fe and Si are interestingly comparable to those of gas-phase abundances in 
intermediate-velocity interstellar clouds which are recently measured using Goddard High 
Resolution Spectrograph, suggesting that resilient cores of grains are not easily destroyed in 
shocks ( |Fitzpatrick 1996| ), We will discuss the implications for grain destruction more in §6.2 
below. 
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Fig. 10. — Spectrum of the P II 60.6 \im line toward 3c 391:BML. The solid curve is a Gaussian 
fit to the line with width fixed to the LWS grating spectral resolution of AA = 0.29 \im. 

Phosphorus — The abundance of P is known for only a few lines of sight in the interstellar 
medium. Our abundance estimate in the post-shock gas uses the observed brightness of the 60.6 
fim line, for which this is the first reported detection (at the time of writing). Figure shows the 
spectrum of this newly-detected line. We used the radiative decay rates are as given by Mendoza 
& Zeippen (1982), and we use the collision strength 0( 3 Pi — ► 0) = 1.46 flKrueger fc Zsyzak 197C| ). 
The P abundance relative to O is very well constrained because the energy levels of the two 
observed ions are similar, and both elements were observed in their dominant ionization state. We 
find that the observed line brightness is consistent with solar abundance of P. Along the line of 
sight toward ( Oph, it appears that about half of interstellar P is normally locked in grains 
( Savage fc Sembach 1996] ). Our observation suggests that some or all of the P-bearing solid 
material was vaporized by the strong shock wave, as we find for Fe and Si (see §|6.2j). 



5. Shock-excited molecules 
5.1. Ho 



The S(3) and S(9) lines of H2 were both detected from each of the supernova remnants we 
observed, and the line brightnesses are listed in Table ^. These H2 lines arise from shocks in the 



molecular gas, like the CO, OH, and H2O lines (Reach & Rho 1998). Using our cartoon model, 



we can now localize which region(s) give rise to which lines and further constrain the properties of 
the dense clumps (region C). First, for region M, if we assume that the H2 rotational lines are 
collisionally excited with rotational excitation temperature equal to the kinetic temperature, we 
predict that the brightest H2 line is the S(3) line. The brightness of the S(3) line from region M, 
assuming all of the H is in H2, is 2 x 10 -4 erg cm -2 s _1 sr _1 . This is in very good agreement with 
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the observed line brightness, suggesting that the H is largely molecular behind the region M 
shocks. The column density of H2 is just enough that H2 molecules are self-shielded from 
dissociating photons in the interstellar radiation field ( |Hollenbach, Werner, fc Salpeter 1971 ). 



We cannot explain the brightness of the higher H2 S(9) rotational line, or the higher- lying 
ro-vibrational lines in the near-infrared, using the region M shocks. These lines most likely arise 
in molecular regions that are much warmer than the recombined molecules behind a fast, 
dissociative shock. Instead, the S(3) lines probably arise from non-dissociative shocks into denser 
gas (discussed below). 



In a previous paper ( Reach &, Rho 19~9g| ), we derived the physical properties of region C using the 



excitation of CO, OH, and H2O. In that paper, we found the puzzling result that abundance of 
H2O was lower than expected in a region with 'warm' chemistry, where O is efficiently converted 
into molecules and OH is converted into H2O. Two papers have shown that OH can be formed 
from dissociation of H2O ( Lockctt et al. 1998 , Wardlc 1999j ), so the abundance OH is no longer 



such a puzzle. We now reexamine the H2O abundance relative to H2. The column density of H2O 
was found to be ~ 3 x 10 17 cm~ 2 using the far-infrared line brightnesses. The column density of 
H2 was taken from the volume density times path length, yielding ~ 8 x 10 23 cm~ 2 . We suggested 
this column density was consistent with the brightness of the H2 S(3) line, but we now suspect 
that much of the S(3) line arises from region M and not the dense shock of region C. If this is 
true, then the H2 column density of region C could be much lower than we estimated before. Let 
us start with some assumptions about the molecular abundances and infer the H2 abundance. 
First, CO:OH:H20 is 100:1:15 from our observations ([Reach fc Rho 199~8| ). Assume that the C 
and O are largely in molecules. If a fraction fco of the C is locked in CO, then the abundance 
ratio of H2O/H2 is Xofco, where Xo is the abundance of O. The column density of H2 is then 
Ar(H 2 )( c ) = 3 x 10 21 /co- Such a column density for region C seems reasonable, but it is lower 
than our previous estimate based on the observed size of the CO emitting region in the millimeter 
observations ( Reach fc Rho 1999| ). This could be explained by a low filling factor for the dense 



gas, of order 0.007/^q in the ISO SWS beam. Taking fco = 0.5 as an example, the size of the 
emitting region would be of order 1", and it would be highly structured in future, 
higher-angular-resolution observations. If true, this new interpretation would change our estimate 
of the mass of the 3C 391:BML post-shock gas, making the self-gravity of the shocked gas 
negligible, and weakening our argument that the dense clump might imminently collapse to form 
stars (Reach & Rho 1999). Table |5| shows the cartoon model for region C assuming fco = 0.5. 



The H2 0-0 S(9) line cannot be assigned to any of the emitting regions that we have defined so 
far. It's energy level is too high above the ground state to be excited in the cooling region where 
the far-infrared lines are produced. It is unlikely that the S(9) line can be produced in the region 
M shock, even upstream from the cooling region, because the H2 is probably still dissociated. 
Thus we suspect that the S(9) line is produced upstream from the region C shock. To get the 
observed line S(9) brightness, in an emitting region with rotational excitation temperature of 
order 1500 K, a column density of order 10 19 cm -2 is needed, which could easily fit in the 
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upstream region considering the cooling region has a column density 600 times larger. The 1500 



K excitation temperature seems high, but high excitation is predicted by shock models (Draine et 



al. 1983| ) and is required to explain observed line ratios in IC 443 ( |Rho et al. 2000| , Richter 



Graham, fc Wright 1995] ), Future observations should be able to distinguish which shocks produce 
the S(3) and S(9) lines, by measuring the width of the lines. According to our cartoon model, 
much of the S(3) line comes from region M and would have a width of 100 km s _1 . The S(9) line 
comes upstream from region C and would have a width of 30 km s^ 1 . The line profiles of 
progressively higher excitation lines would evolve from wide to narrow lines. 




Fig. 11. — Spectrum of an unidentified line toward 3C 391:BML. The solid curve is a Gaussian fit 
to the line with width fixed to the LWS grating spectral resolution of AX = 0.29 \xm. 



5.2. Unidentified line at 74.3 /zm 



One relatively bright line remains unidentified. Figure 11 shows the spectrum of the 74.3 /jm line 
for 3C 391:BML. The line was present in our complete LWS spectra of W 28, W 44, and 3C 391, 
and it has also been seen in the remnant RCW 103 ( Oliva et al. 1999b| ) and the planetary nebula 



NGC 7027 ( |Liu et al. 1996] ). Table || lists the brightnesses and central line wavelengths for our 
remnants. The wavelengths are all consistent with 74.26 (im. All of the spectra that contain this 
line contain also contain both atomic fine structure and molecular lines. Based on our list of 
plausible lines compiled above in Tables [3| and |], and a careful inspection of all atomic fine 
structure lines, we can find no candidate atomic or ionic source for the line. (The only fine 
structure line that matches the observed wavelength is a transition among high energy levels of 
[Ti III], which is extremely unlikely based both on abundance and energetics.) Therefore, we 
suspect this line is molecular in origin. Despite searching the JPL spectral line database ( [Pickett 
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et al. 1996), we cannot identify the origin of this line, which is apparently caused by an unusual 
molecule or meta-stable state not yet found in laboratories or planetary atmospheres. 

6. Continuum and Grain Destruction 

6.1. Continuum spectrum 

The moderate-resolution far-infrared spectra are dominated by a bright continuum underlying the 
spectral lines. Because the remnants we discuss here are in the galactic plane, much of this 
continuum emission is likely to be unrelated to the remnant, produced instead by clouds along the 
line of sight both in front of and behind the remnant, as well as the molecular cloud close to the 
remnant. We have attempted to extract the remnant contribution to the continuum emission by 
searching for distinct spatial and spectral variations. In this paper we present a spectral 
separation of the remnant emission from the unrelated emission; in paper 1 and in future work, 
we use the reference positions and spatial correlations with the O I line brightness^]. 



3 We would like to clarify an issue from paper 1. Despite a statement by Oliva et al. (1999a) that our continuum for 
W 44 is 10 times too bright compared to IRAS data, the continuum brightness of the remnant as observed with ISO 
appears real. The IRAS data suffer from lower angular resolution, which makes it extremely difficult to separate the 
remnant from unrelated galactic emission. The continuum emission will be discussed in more detail in a forthcoming 
paper. 



Table 8: Unidentified line 



source 


wavelength 


Line Brightness 




(jmx) 


(10~ 4 erg s _1 cm _2 sr _1 ) 


3C 39LBML 


74.26 ± 0.03 


0.31 ±0.05 


W 44 


74.32 ± 0.06 


0.11 ±0.03 


W 28 


74.25 ± 0.03 


0.26 ±0.04 
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Fig. 12. — Far-infrared continuum spectra toward molecular shocks in (a) W 28, (b) W 44> an d ( c ) 
3C 391. The jagged curve shows the observed brightness. The solid curve shows the spectral model 
for unrelated dust along the line of sight. This model consists of two components of dust associated 
with atomic gas (T^ ~ 20 K) and molecular gas T^ ~ 13 K); each component is shown as a dotted 
curve. The hatched region shows the excess emission that could be the continuum emission from 
the remnants. The two dotted curves show the individual contributions of unrelated atomic and 
molecular gas. 
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The continuum spectrum of 3C 391:BML is shown in Figure [T^. Even by inspection, it is evident 
that the spectrum cannot be fitted with a modified blackbody: there is a significant positive 
excess at wavelengths shorter than 120 /im. It is not surprising for the continuum spectrum to 
require multiple components. In the COBE FIRAS observations of the galactic plane in a 7° 
beam, the spectrum of the inner galaxy at wavelengths of 100-200 /im has contributions from two 
components: 20 K dust associated with atomic gas, and 13 K dust associated with molecular gas 
( Roach et al. 1995] , Lagache et al. 199S ). At wavelengths shorter than 100 /im, continuing all the 
way to the PAH features in the mid-infrared, emission from small, transiently-heated dust grains 
becomes important ( Desert, Boulanger, fc Puget 1990| ) . For our supernova remnant observations 
in an 80" beam, we expect all components to be brighter than in the COBE/FIRAS data, because 
(1) the atomic gas column density peaks sharply in the galactic plane, (2) the lines of sight pass 
directly through giant molecular clouds, and (3) the emission from the remnants themselves are 
completely beam-diluted in the COBE FIRAS beam. 

To estimate the remnant contribution to the spectrum, we used a nominal spectrum for dust in 
the local interstellar medium, which big grains, at a temperature of 17.5 K, and small grains with 
a brightness that matches the COBE Diffuse Infrared Background Experiment ( Reach fc 



Boulanger 1998| and references therein). We use two components, for the atomic and molecular 



gas, with radiation fields equal to x times the local interstellar radiation field. For each 
component, the temperature of big grains scales as and the brightness of the small grain 
emission scales as x- The atomic regions were assumed to have x = 2 (to account for the 
enhanced radiation field in the inner Galaxy), and the molecular regions were assumed to have 
X = 0.5 (to account for extinction of the ultraviolet radiation field inside molecular clouds). The 
scale of the model and the ratio of the two components was set to match the far-infrared 
spectrum at wavelengths longward of 120 /im. 



The continuum spectral models are shown for each remnant in Fig. 12. The models with a single 
diffuse ISM component cannot explain the brightness of the continuum at wavelengths shortward 
of 110 /im; there is a clear excess at 40-110 /im. The amount of excess emission depends rather 
sensitively on the assumed spectrum of the unrelated interstellar material, so the present results 
must be used with caution. There is a good correlation of the 74 fim continuum with the O I line 
brightness (paper 1), supporting our claim that a significant portion of the 40-110 /im emission is 
from the remnants. There is also a tentative detection of 3C 391 at 60 /im using IRAS data 
([Arendt 19961 ). This excess emission has a color temperature of ~ 40 K, corresponding to an 
enhancement of the radiation field by more than 2 orders of magnitude. Such a high radiation 
does not occur in the diffuse interstellar medium over many-pc scales such as observed here. 
Instead, we must be seeing dust that is either physically different from that in the diffuse 
interstellar medium, or dust that is excited locally, for example behind the shock fronts. The 
optical depth of the remnant emission at 100 /jm is r ~ 2 x 10 -4 , and the total surface brightness 
of the far-infrared continuum, assuming a u 2 B^(A2K) spectrum, is 2 x 10 -2 erg s _1 cm" 2 sr" 1 . 
Thus the infrared continuum is significantly brighter than any of the individual spectral lines. 
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The luminosity of the excess emission, within the single LWS beam that we observed for each 
remnant, ranges from ~ 6OOL for 3C 391:BML to ~ WL & for W 44 and W 28; however, these 
numbers are rather uncertain. Using the optical depth required to produce the excess emission, 
the opacity per unit mass for 'normal' interstellar dust (Draine & Lee 1984), the mass of 
shock-excited dust within the single LWS beam that we observed for 3C 391:BML, is ~ LM©. 



6.2. Grain destruction 

In order for the Fe and Si emission lines to be as bright as observed, at least part of the refractory 

dust mass must have been vaporized. In the cool diffuse cloud toward £ Oph, 95% of Si and 

> 99% of Fe are locked in dust grains ( Savage fe Sembach 1996| ). Using our single-slab model, we 



found that 50% of the Si and 30% of the Fe are required to be in the post-shock gas, which would 
mean that about half of the dust mass is destroyed by the shock. The single-slab model will, 
however, give less accurate measures of the abundances for Fe and Si relative to O (for which the 
model parameters were determined), because the Fe and Si lines we observed are from higher 
energy levels than those of O. There is likely to be a stratified post-shock emitting region, with 
somewhat warmer regions producing some of the Si and Fe lines. Thus, the exact fraction of 
refractory dust mass that was vaporized is highly uncertain. We have experimented with a 
multiple-slab model, using the approximate density and temperature profiles behind a 100 km s _1 
J-shock into 10 3 cm -3 gas ( |Hollenbach fe McKee 1989| ), to see how much of the observed 35.2 and 
26 /um line brightnesses could be produced by warmer layers closer to the shock. The multiple 
slab model predicts relatively brighter 26 fim lines for a given column density, by up to a factor of 
two in extreme cases. Thus it is possible that the gas-phase abundance of Fe was overestimated 
by up to a factor of 2 using the single-slab model, meaning that as little as 15% of the Fe is 
required to produce the observed lines. An identical result is obtained for Si. Even with the lower 
gas-phase abundance of Si and Fe, we can conclude that a non-trivial fraction of the Fe-bearing 
dust mass was vaporized. It is essentially impossible to produce all of the observed brightness of 
the 26 /im line from normally-depleted pre-shock interstellar matter. 

Using the observed traces of the vaporized and surviving dust grains, we can work backwards to 
estimate the fraction of dust that was destroyed. The total mass of Si and Fe vapor in 
3C 39LBML is about 0.11 M Q , within the LWS beam. Based on the dust-phase elemental 
abundances in a cold interstellar cloud, the total mass of elements (mainly O, Fe, C, Si, and Mg) 
in dust is 3.6-5.0 times the mass of Si and Fe in dust, with the lower number coming from B-star 
reference abundances and the higher number coming from Solar reference abundances ( [Savage & 



Sembach 1996j ). Therefore, we infer that about 0.5 M@ of dust was vaporized by shocks into 
3C 39LBML. For comparison, our far-infrared conntinuum observations indicated at least 1 M & 
of dust remaining in solid form after the shock. Therefore, we estimate that approximately 1/3 of 
the dust mass was vaporized by the shocks. 



How does our inferred fraction of dust vaporization compare to theoretical predictions? Jones, 
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Tielens, & Hollenbach (1996) have calculated the destruction probabilities for dust grains in shock 
fronts with a range of shock velocities and pre-shock densities. Their models included sputtering 
and shattering: sputtering removes material from the grain via collisions with the hot gas, while 
shattering destroys large grains and creates small grains from them. The models for v s = 100 
km s _1 are most relevant for the shocks we discuss in this paper. Clearly there are faster shocks 
that create the X-ray-emitting plasma, and slower shocks that create the molecular line emission. 
But the ionic lines from which we determined the abundances arise from a layer with a post-shock 
density of order 10 3 cm -3 (§||) into which shocks with a velocity of about 100 km s _1 have been 
driven (§ |2.3j ). The model of Jones et al. (1996) with parameters most similar to our present case 
has pre-shock density tiq = 25 cm -3 . While their models do not extend to higher no, the 
dependence on no is surprisingly mild, with a slight trend toward more dust destruction for 
denser shocks. The models of Jones et al. (1996) predict that about 37% of silicate dust and 12% 
of graphite dust are destroyed, mostly via non-thermal sputtering. These predictions are in 
excellent agreement with our results. 



7. Nature of the shock fronts powering the emission lines 

There are at least two significant theoretical models that predict the structure (density and 
temperature) and emission behind radiative shocks comparable to those that seem to be occurring 
in the molecular supernova remnants we observed. One model, by Hollenbach & McKee (1989; 
hereafter HM89) , includes an extensive network of chemical reactions and detailed treatment of 
the dynamics of the post-shock gas (from [Hollenbach McKee 1979 ). The other model, by 



Hartigan, Raymond, & Hartmann (1987; hereafter HRH87), includes more ionization states and 
more ultraviolet and optical lines. We will use the HM89 model as a benchmark, and we will use 
the [O I] 63 fim line as the normalizing factor. In the HM89 model, the 63 line is a significant 
coolant and should reliably indicate the total energy passing through the shock. Based on the 
observed width of the line (§^^), we consider the models with shock velocity V s = 100 km s _1 . 
The HM89 model at that velocity that matches the observed [O I] 63 fim brightness has pre-shock 
density no = 10 3 cm~ 3 . Such a pre-shock density is reasonable considering that we are observing 
shocks into molecular clouds. This shock corresponds to what we called 'region M' (§4.2), 



although we suggested a somewhat lower pre-shock density of 10 2 cm -3 for this shock. In fact we 
do not directly measure no and instead estimated it from the density in the emitting region, so a 
higher no could be consistent with the observations. The compression by the shock depends on 
many unknown factors, including the strength of the preshock magnetic field. 

Before stepping through the comparison of more observed results to the models, we should clarify 
that the models are themselves significantly different. The physics of shock fronts, especially into 
dense gas and dust clouds, with so many modes of energy release and interaction, is exceedingly 
complex and cannot be expected to be accurately predicted by a simple theoretical model. Let us 
compare the HM89 and HRH87 models for n = 10 3 cm" 3 and V s = 100 km s _1 . The HRH87 
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model is called D100 in their paper. The magnetic field is negligible in the D100 model, but 
comparing the E100 and B100 models it appears that a stronger magnetic field makes relatively 
little difference in predicted line brightness within the context of the HRH87 model. Comparing 
the line brightnesses between the HM89 and HRH87 models, it is clear that there are very 
significant differences. An extreme case is the [O I] 6300 A line, which is very bright in the HM89 
model but weak in the HRH87 model; the difference is a factor of about 8. A significant part of 
the difference could be neglect of the higher ionization states in the HM89 model; we suspect the 
[O II] and [O III] optical lines are much brighter than [O I]. Other significant differences are the 
importance of Lya, many high-excitation ultraviolet lines, and the 2-photon continuum; all of 
these are very bright in the HRH87 model. Part of the difference could be neglect of the radiative 
transfer through dust in the HRH87 model; we expect most of the ultraviolet photons are 
absorbed by the cooler layers of the post-shock gas and dust or by the remaining unshocked gas 
and dust. 

Comparing our observations to the theoretical shock models, we reach the following conclusions: 

Dust emission — We observe a far-infrared continuum toward the shock fronts, in excess of what 
we would expect from dust in the quiescent pre-shock gas. The dust has a color temperature 
distinctly warmer than dust in molecular or atomic clouds. This could be due to significant 
collisional or radiative heating in the post-shock gas raising the equilibrium temperature of the 
grains, or it could be due to a significant difference in the size distribution of the post-shock 
grains. Indeed the post-shock size distribution is predicted to be significantly different because 
the larger grains are shattered (|Jones, Tielens, fe Hollenbach 19"96 ). The HM89 model predicts 



too little [Fe II] 26 fim and [Si II] 35 /jm emission by about a factor of 10. This is probably due to 
treatment of dust destruction: the elemental abundances in the HM89 model are set at depleted 
cosmic abundances typical of what may occur in the pre-shock gas, while dust destruction is now 
predicted to be significant at these shock velocities ( pones, Tielens, fe Hollenbach jggg ). The 
HM89 model predicts some infrared continuum, but its flux is weaker than the [O I] 63 /mi line, 
which does not agree with our observations: we find that the continuum flux is much greater than 
that of any of the emission lines. The HRH87 is on the opposite extreme: they use the entire 
cosmic abundance in the gas phase, implicitly assuming complete destruction of the dust and 
predicting no infrared continuum. 

Infrared versus ultraviolet — The infrared emission escapes from the shocks and arrives at our 
telescopes, while the ultraviolet (and optical) emission either never escapes the shock or is 
absorbed before reaching us. The HRH87 model predicts copious ultraviolet emission, with the 
strongest emission (in order of brightness) arising from Lya, 2-photon continuum, C III 977 A, and 
He II 304 A. The sum of the 5 infrared lines they report is still 20 times weaker than Lya and 3 
times weaker than just the He II 304 A line. Thus ultraviolet lines (and 2-photon continuum) are 
by far the dominant cooling behind radiative shocks according to the HRH87 model. On the other 
hand, the HM89 model predicts a significantly different spectrum, with relatively brighter optical 
and infrared lines. The brightest line in the HM89 model is the optical [O I] 6300 A line, followed 
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by Ha and comparable brightnesses of infrared [O I] 63 /xm, optical [S II] 6731 A and [N II] 
6560 A, and ultraviolet C II] 2326 A lines. Remarkably, the brightness of the C II] 2326 A and 
[O II] 3726+3729 A lines, and Ha and H/3, agree for the two models, but there are many brighter 
ultraviolet lines in the HRH87 model. The difference between the two models is at least partially 
due to the different assumptions about dust grains. The HRH87 model has no dust grains, so that 
ultraviolet transitions that are optically thick can scatter from ion to ion and eventually escape 
the shock or have a quantum decay into two photons. If dust grains were included, the ultraviolet 
photons would not last long before being absorbed. The ultraviolet photons are produced in a 
thin layer behind the shock, and their energy would heat the surrounding g radiative 
precursor into the pre-shock gas and as an enhanced radiation field in the cooler layers of the 
post-shock gas. In the HM89 model, the radiative transfer does take into account absorption of 



the ultraviolet lines by dust grains (Hollenbach & McKee 1979), which explains why this model 



predicts much weaker ultraviolet emission than the HRH87 model. Our detection of a far-infrared 
continuum suggests that the models must include the effects of dust on the ultraviolet photons. 



Molecules — We detected both the S(3) and S(9) lines of H2 (§ f2.2p , as well as OH and weak H2O 
emission ( Reach &: Rho 1999 ). The model of HRH87 does not include molecules, so we will not 



discuss it here. The model of HM89 includes a chemical reaction network, and they predict the 
brightness of the S(3) and S(9) lines to be 9 x 1(T 6 and 8 x 10 -5 erg s 1 cm 2 sr , respectively. 
The predicted brightness of the S(9) line is about half of what we observed, suggesting that a 
significant fraction of the S(9) line could be produced by a J-type shock. However, the predicted 
brightness of the S(3) line is far lower than we observed. Indeed, for all 4 of the molecular shocks 
we observed, the S(3) line is brighter than the S(9) line; this is contrary to the model, which 
predicts S(3) about 10 times weaker than S(9). Instead, it appears that a denser, non-dissociating 
shock is required to explain the S(9) line brightness. This picture agrees with presence of a 
preshock cloud with a range of densities, as shown in our cartoon model. Therefore, it appears 
that the shock models cannot simultaneously predict both the ionic lines and (all of) the 
molecular lines with no = 10 3 cm -3 and V s = 100 km s . Instead, it appears that a denser, 
non-dissociating shock is required to explain the S(3) line brightness. 

To explain the brightness and excitation of the OH and H2O emission, and the ~ 30 km s _1 
width of the CO and CS lines, we suggested C-type shocks into gas with density > 10 4 cm -3 
(Reach &: Rho 98, 99). Models for such shocks were developed by Draine, Roberge, & Dalgarno 
(1983; hereafter DRD). For n = 10 4 cm" 3 and V s = 30 km s -1 , the DRD models predict that H 2 
is the dominant coolant, with the brightest lines being S(3) and S(5). In the DRD model, the S(9) 
line is about 10 times weaker than the S(3) line. Given that we observe S(3) and S(9) lines of 
comparable brightness, it appears that we are required to have the bulk of the the S(3) line arise 
from a C-type shock and part of the S(9) line arise from a J-type shock. The observed surface 
brightness of the S(3) line is lower than the predicted brightness from the DRD model 
(Ipred = 3 x 10 -3 ), but that could be due to beam dilution of small, dense clumps that produce 
the C-type shock. 



-36- 



Comparison to a model of W 44 — Two recent papers by Cox et al. (1999) and Shelton et al. 
(1999) presented a model for the remnant W 44. In this model, all observed properties of the 
remnant, from radio emission to 7-rays, and specifically including the infrared lines, are produced 
by a blast wave into a smooth medium with a density of 6 cm -3 . This model is in marked 
contrast to the one we presented here, with pre-shock densities ranging from < 1 to > 10 2 cm -3 , 
and some discussion may help to resolve confusion about how they could both be explaining the 
same data. First, we note that the evidence for molecular cloud interaction is very strong. 
Shocked CO line with widths of ~ 30 km s _1 were revealed by Seta et al. (1998) and ourselves 
(this paper, Fig. [?]). Second, we note that the uniform pre-shock conditions in the Cox et al. 
model are not realistic in a large region of the multiple-phase interstellar medium. 

In our model, we attribute the molecular line observations of W 44, including H2 near-infrared 
emission ( [Beta et al. 1998 ), infrared CO lines ( [Reach fc Rho 1998| ), wide millimeter CO and other 



lines ( Seta et al. 1998] ), and radio OH 1720 MHz masers ( plaussen et al. 1997] ), to dense, shocked 



gas. These lines cannot be produced by an ionized low-density post-shock region. The analogy to 
IC 443 is useful — W 44 is very similar in many regards, with similar brightnesses of the 63 /im line 
and brightnesses and widths of the millimeter-wave CO (and other) lines — and few would deny 
that IC 443 is interacting with a molecular cloud. Finally, the observed brightness of the [O I] 63 
fim line cannot easily be explained by a low-density medium, in contrast to the discussion claimed 
by Cox et al. (1999). The excitation of [O I] evidenced by the 145/63 /im ratio requires a higher 
density region (Fig. [8]). Obviously, [O I] also requires a neutral region, for if the H is ionized the O 
will also be ionized. The observed lines were so 'shockingly' bright that we concluded immediately 
that they were due to a dense shock ( Reach Sz Rho 1996| ). In this paper we considered a wide 



range of densities and ruled out the low densities because (1) they are inconsistent with the line 
ratios, and (2) the path length required through a low-density gas would be longer than the 
remnant size in order to build up enough column density to produce the observed line brightness. 
Cox et al. reduced the brightness of the 63 [im by a factor of 100 to match their models. One of 
the reduction factors is a a factor of 10 to match the estimated average over the entire remnant. 
Another factor applied by Cox et al. was to enhance their prediction because of limb brightening 
and local rippling of the shell. While this is perfectly reasonable, this correction was applied to 
the low-resolution 63 //m line observations, which trace relatively cool gas, and not to the 



higher-resolution Ha observations ( Giacani et al. 1997 ), which trace warmer gas that might be 



even more edge brightened. All of these factors probably cancel out to first order when the 
average over our large beam is made, so that the [O I] /Ha ratio is not matched by their model. 

The truth is probably that there were both high and low density regions in the interstellar 



medium around the W 44 progenitor star (as considered in this paper and by Chevalier 1999). 
The particular lines of sight studied in this paper are centered on a special regions, where 
evidence for a dense shocks are overwhelming. Away from these positions, it is likely that the 
lower-density shocks contribute some of the low-level 63 /im emission. In Figure ||, we are pointed 
right at a region C clump, while the average over a large region would include different regions in 
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different proportions. Because the 63 /im line brightness is expected to scale approximately as 
uqVs (HM89), the denser regions contribute relatively more to the observed brightness, when they 
are present. Further observational work may be able to separate the dense shocked gas from the 
more rarefied and widespread interclump gas, using imaging to resolve the clumpy spatial 
distribution and spectroscopy to separate the narrower line widths. 



8. Conclusions 

Based on 750 spectroscopic observations of supernova remnants interacting with molecular 
clouds, we found evidence for a wide range of pre- and post-shock conditions; these properties are 
summarized in Table [| and illustrated in Figure ^. Of the three density regimes that we identified, 
most of the infrared emission arose from the shocks into moderate-density molecular gas ('region 
M'). The energy of these shocks is radiated via continuum emission from surviving dust grains 
and [O I] 63 /mi and [Si II] and other infrared atomic fine structure lines. The dust continuum 
was difficult to separate from the far-infrared emission from cold, unrelated dust, but the shocked 
dust was evident at 80-100 /xm as a warmer component of the spectrum. Some of the dust grains 
must have been destroyed, because we observe bright emission lines from dust 'vapors,' including 
the [Fe II] 26 /im line. To produce the observed line dust vapor lines, we require that 1/3 of the Si 
and Fe-bearing dust mass was vaporized in the shocks. Using the Fabry-Perot observations, the 
width of the 63 /jm line was found to be ~ 100 km s . Theoretical models can reproduce the 
brightness of the 63 /xm line with such shock velocities and the properties of 'region M' 



( Hollenbach fe McKee 1989 ); theoretical models can also reproduce the observed amount of grain 



destruction ( Jones, Tielens, fe Hollenbach 1996 ). Higher-density gas was required to explain the 



bright H2 line emission, leading to our 'region C,' for clumps. Lower-density gas was required to 
explain the higher-ionization lines, leading to our 'region A,' for atomic gas. 

A general conclusion from these observations is that molecular shock fronts are copious producers 
of infrared emission. We anticipate many applications of these types of observations for further 
understanding the nature of molecular supernova remnants, such as the fate of dense clumps, the 
fate of dust grains, and the effect of cooling on the remnant evolution. In turn, these studies can 
elucidate, to some extent, the nature of the molecular clouds before the shocks. If molecular 
clouds were uniform, or nearly so, then we would expect an orderly progression of pre-shock to 
post-shock gas, with radiative coolants characteristic of a single shock velocity. On the contrary, 
we observe emission from coolants as diverse as multi-atom molecules, multiply-ionized atoms, 
and dust grains and their vapors. 



We would like to thank Pierre-Olivier Lagage for helpful discussions in interpreting atomic fine 
structure lines, and we thank Emmanuel Caux, Steve Lord, and Sergio Molinari for clarifying 
details of the LWS instrument and calibration. 
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